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ABSTRACT 

We present a study of the redshift evolution of the projected correlation function of 593 X-ray 
selected AGN with Iab < 23 and spectroscopic redshifts z < 4, extracted from the 0.5-2 keV X-ray 
mosaic of the 2.13deg 2 XMM-COSMOS survey. We introduce a method to estimate the average bias of 
the AGN sample and the mass of AGN hosting halos, solving the sample variance using the halo model 
and taking into account the growth of the structure over time. We find evidence of a redshift evolution 
of the bias factor for the total population of XMM-COSMOS AGN from b{z = 0.92) = 2.30 ± 0.11 to 
b(z = 1.94) = 4.37±0.27 with an average mass of the hosting DM halos logM^n^M®] ~ 13.12±0.12 
that remains constant at all z < 2. 

Splitting our sample into broad optical lines AGN (BL), AGN without broad optical lines (NL) and 
X-ray unobscured and obscured AGN, we observe an increase of the bias with redshift in the range z = 
0.7 — 2.25 and J — 0.6 — 1.5 which corresponds to a constant halo mass logM^n^M^] ~ 13.28 ± 0.07 
and logMolh^M®} ~ 13.00 ± 0.06 for BL /X-ray unobscured AGN and NL/X-ray obscured AGN, 
respectively. 

The theoretical models which assume a quasar phase triggered by major mergers can not reproduce 
the high bias factors and DM halo masses found for X-ray selected BL AGN with Lbol ~ 2 x 10 45 erg 
s _1 . Our work extends up to z ~ 2.2 the z < 1 statement that, for moderate luminosity X-ray selected 
BL AGN, the contribution from major mergers is outnumbered by other processes, possibly secular 
such as tidal disruptions or disk instabilities. 

Subject headings: Surveys - Galaxies: active - X-rays: general - Cosmology: Large-scale structure of 
Universe - Dark Matter 



1. INTRODUCTION 

Investigating the clustering properties of active galac- 
tic nuclei (AGN) is important to put tight constraints 
on how the AGN are triggered and fueled, to identify 
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the properties of the AGN host galaxies, and to under- 
stand how galaxies and AGN co-evolve. In addition, in 
the framework of the cold dark matter (CDM) structure 
formation scenario, clustering properties or the bias of 
AGN, may be related to the typical mass of dark mat 
ter (DM) halos in which they reside Mo fc White) 1 199 
Shet h fc Tormenl 119991 : iSheth et all 120011: iTinker et al 



2005) and allow various types of AGN to be placed in 



a cosmological context 

Recently, several studies have been made, employ- 
ing spectroscopic redshifts to measure the three dimen- 
sional correlation function of X-ray AGN. The major- 
ity of the X-ray surveys agree with a picture where X- 
ray AGN are typically hosted in DM halos with mass 
of the order of 12.5 < logM DM \hr Y M q] < 13.5, at low 
(z < 0.4) and high (z ~ 1) redshi ft dGilli et all [2005 : 
Yang et all [20061: iGilli et~aL1 [200l iHickox et all 120091: 



Coil et all 120091; IKrumpe et all 120101: ICappelluti et al l 
20101 ). This implies that X-ray AGN more likely reside 



in massive DM halos and preferentially inhabit dense en- 
vironment typical of galaxy groups. 

There have been attempts to detect X-ray luminos- 
ity dependence of t he clustering. A t z ~ 1, neither 
IGilli et al.l (|2009l ) nor lCoil et aT1(|2009l) found significant 
dependence of the clustering amplitudes on the optical 
luminosity, X-ray luminosity or hardness ratio, partially 
due to the lar ger st atistical errors. Recent works by 
IKrumpe et all |2010l) and ICappelluti et al.l (|2010l ) found, 
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however, that high X-ray luminosity AGN cluster more 
strongly than low X-ray luminosity ones at 2a level for 
z ~ 0.3 and z ~ 0, respectively. 
Until recently, the clustering of AGN has been studied 



as 2d F (2QZ, iCroom et all 



mainly in optic al, particularly in la r ge area surveys such 

" I2005T : IPorci 

(S 

2006T iShen eTa l. 2009; [Ross etaTI [2009fl . ICroom et al ' 



Po rciani fc Norbera 



2006) and Sloan Digi t al Sky Survey ( S PSS. iLi et al 



(2005) analysed the clustering of 2QZ QSO as a func- 
tion of redshift finding a strong evolution of QSO bias, 
with = 0.53) = 1.13 ± 0.18 at low redshift and 

b Q (z = 2.4 8) = 4.24 ± 0.53 at high redshift, as also 
observed in iPorciani fc Norbergj (|2006t ). The evidence 
of an evolution over time of the bias facto r for SDSS 
quasars has been found in IShen et al.l ([2009), with bias 
values ranging from 6q(z = 0.50) = 1.32 ± 0.17 to 
b Q (z = 3.17) = 7.76 ± 1.44. The results from these sur- 
veys have also shown that the bias evolution of optically 
selected quasars is consistent with an approximately con- 
stant mass at all redshifts of the hosting DM halo in the 
range logM DM ~ 12.5 - 13[/i _1 M ]. 
Besides models of major mergers between gas-rich galax- 
ies appear to naturally produce the bias o f quasars as 
a function of L and z (Hopkins et al.1 120081: IShen! 120091: 
Shankar et ail 120091 120101: iShankarl 120101: iBonoli et al.l 
2009f ). supporting the observations that bright quasars 
host galaxies present a preference for merging systems. 
It is still to be verified if the results from optical surveys 
can be extended to the whole AGN population and in 
particular to the X-ray selected AGN. 

In this paper, we concentrate on the study of the bias 
evolution with redshift using different X-ray AGN sam- 
ples and we focus on the estimation of the bias factor 
and the hosting halo mass using a new method which 
properly account for the sample variance and the strong 
evolution of the bias with the time. 
The paper is organized as follows. In section Sj2]we de- 
scribe the XMM-COSMOS AGN sample and the AGN 
subsets used to estimate the correlation function. In $3] 
we describe the random catalog generated to reproduce 
the properties of the data sample and the method to 
measure two-point statistic is explained in f}4] The re- 
sults of the AGN auto-correlation based on the standard 
method of the power-law fitting of the signal and using 
the two-halo term are given in fj5j In fj6] we present our 
own method to estimate the AGN bias factor and the 
DM halos masses in which AGN reside, solving the sam- 
ple variance and the bias evolution with redshift and in 
iJ7]the results. In fJ5]we present the redshift evolution of 
the bias factor and the corresponding DM halo masses 
for the different AGN subsets. We discuss the results in 
the context of previous studies in £[9] and we conclude in 
CTUl Throughout the paper, all distances are measured 
in comoving coordinates and are given in units of Mpc 
h" 1 , where h = H /100 km/s. We use a ACDM cosmol- 
ogy with n M = 0.3, tl A = 0.7, n b = 0.045, a$ = 0.8. 
The symbol log signifies a base-10 logarithm. 

2. AGN CATALOG 

The Cosmic Evolution Survey (COSMOS) is a multi- 
wavelength observational project over 1.4 x 1.4 deg 2 of 
equatorial field centred at (RA, DEC) J2000 = (150.1083, 
2.210), aimed to study AGN, galaxies, large scale struc- 
ture of the Universe and their co-evolution. The sur- 



vey uses multi waveleng th imaging from X- ray to radio 
bands, including HST (IScoville et al.l 120071) . SUBARU 
(jTaniguchi et alj2007Q, Spitzer (jSanders et alj2007h and 
GALEX (jZamoiski et all |2007[) . The central 0.9 deg 2 
of the COSMOS field has bee n observed in X-r ay with 
Chandra for a total of 1.8 Ms (| Elvis et al.ll2009fl . In ad- 
diction spectroscopic campaigns have been carried out 
with VIMOS/VLT and extensive spectroscopic follow- 
up have been granted with the IMACS/Magellan, MMT 
and DEIMOS/Keckll projects. 

XMM-Newton surveyed 2.13 deg 2 of the sky in the COS- 
MOS field i n the 0.5-10 keV energy band for a total of 
~ 1.5 5 Ms (jHasinger et al.l [20071: ICappelluti et all 120071 
2009) providing an unprecedently large sample of point- 
like X-ray sources (1822). 

The XMM-COSMOS catalog has bee n cross-correlated 
with the optical multiband catalog (ICappelluti et al.1 
120071 ). the K-band catalog dMcCracken et al.l 1201 " 



the IRAC catalog (jSanders et all 120071: lllbert et all 



201Ch and the MIPS catalog dLe Floc'h et al.l 12009ft 



Brusa et all (|2010l ) presented the XMM-COSMOS mul- 



tiwavelength catalog of 1797 X-ray sources with opti- 
cal/near infrared identification, multi wavelength prop- 
erties and redshift information (from iLillv et al.l 1 2007 



2009HTrump et al.ll2007tlAdelman-McCarthv et aLl l200€ 
Prescott et alj|2006t iSalvato et al.ll2009D . 



In this paper we focus on the clustering analysis of 
1465 XMM-COSMOS AGN detected in the energy band 
0.5-2 keV, for which we have a spectroscopic complete- 
ness of ~ 53% (780/1465). From this sample of 780 
objects we selected 593 sources with Tab < 23 (this 
magnitude cut increases the spectroscopic completeness 
to about 65%) and redshift z < 4. The redshift dis- 
tribution of the AGN sample (Fig. Q] left panel) shows 
prominent peaks at various redshifts, z ~ 0.12, z ~ 0.36, 
z ~ 0.73, z ~ 0.95, z ~ 1.2, z ~ 2.1. In particular, 
the structure at z ~ 0.3 6 was also obser ved at other 
wavelengt hs in COSMOS (ILillv et"aT1l2007ft and already 
discussed (Gilli et al. 2009). The median redshift of the 
sample is < z >= 1.22. 

The sources have been classified in iBrusa et al.l (|2010f ) 
in broad optical line AGN (BL AGN, 354), non- broad 
optical line AGN (NL AGN, 239) using a combination 
of X-ray and optical criteria, motivated by the fact that 
both obscured and unobscured AGN can be misclassified 
in spectroscopic studies, given that the host galaxy light 
may over shine the nuclear emission. Fig. [2] shows the 
redshift distribution of BL AGN with < z >= 1.55 and 
NL AGN with < z >= 0.74. 

We also studied the clustering properties of X-ray un- 
obscured and obscured AGN derived on the basis of the 
observed X-ray hardness ratio and corrected to take into 
account the redshifts effects. In particular we used the 
hard X-ray band (2-10 keV) (which allows us to sam- 
ple the obscured AGN population) to select a subset of 
184 X-ray unobscured sources (X-unobs hereafter) with 
logNji < 22 cm" 2 and 218 X-ray obscured (X-obs here- 
after) sources with logNn > 22 cm~ 2 , The median 
redshift of the two sub-samples are < z >= 1.12 and 
< z >— 1.30, respectively (see fig. [2j right panel). The 
47% (40%) of BL (NL) AGN have been also observed in 
the hard band and classified as X-unobs (X-obs) AGN. 
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Fig. 1. — Left panel: Redshift distribution of 593 AGN (gold filled histogram) in bins of Az = 0.01, with median z = 1.22. The solid 
black curve is the Gaussian smoothing of the AGN redshift distribution with a z = 0.3, used to generate the random sample (red empty 
histogram). Right panel: distribution of AGN pairs in redshift bins Az = 0.01. 




Fig. 2. — Left Panel: Redshift distribution of XMM-COSMOS AGN (open histogram) selected in the soft band, compared with the redshift 
distribution of BL AGN (blue histogram, upper right quadrant) and NL AGN, (red, upper left quadrant). Lower quadrants show the redshift 
distribution of the random catalogs (open black histograms) for both the AGN sub-samples, obtained using a Gaussian smoothing (gold 
lines) of the redshift distribution of the real samples. Right Panel: Redshift distribution of unobscurcd (dark blue histogram) and obscured 
(magenta histogram) AGN selected in the hard band according with the column density (upper quadrants). Lower quadrants show the 
redshift distribution of the random catalogs (open black histograms) for both the AGN sub-samples, obtained using a Gaussian smoothing 
(gold lines) of the redshift distribution of the real samples. 

3. RANDOM CATALOG 

The measurements of two-point correlation function 
requires the construction of a random catalog with the 
same selection criteria and observational effects as the 
data, to serve as an unclustered distribution to which 
to compare. XMM-Newton observations have varying 
sensitivity over the COSMOS field. In order to create 
an AGN random sample, which takes the inhomogeneity 
of the sensitivity over the field into account, each simu- 
lated source is placed at random position in the sky, with 
flux randomly extracted from the catalog of real sources 



fluxes (we verified that such flux selection produces the 
same results as if extracting the simulated sources from 
a reference input logN-logS). The simulated source is 
kept in the random sample if its flu x is above the sen- 
sitivity map value at that position (Mi vaii et al.l 120071 : 
ICappelluti et al.ir2009f ). Placing these sources at random 
position in the XMM-COSMOS field has the advantage 
of not removing the contribution to the signal due to an- 
gular clustering. On the other hand, this procedure does 
not take into account possible p ositional biases r elated 
to the optical follow-up program. iGilli et al. (2009[), who 
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instead decided to extract the coordinates of the random 
sources from the coordinate ensemble of the read sam- 
ple, showed that there is a difference of only 15% in the 
correlation lengths measured with the two procedures. 

The corresponding redshift for a random object is as- 
signed based on the sm oothed redshift d istribution of 
the AGN sample. As in iGilli et all (|2009| ) we assumed 
a Gaussian smoothing length a z — 0.3. This is a good 
compromise between scales that are either too small, thus 
affected by local density variations or too large and thus 
oversmooth the distribution (our results do not change 
significantly using a z = 0.2 — 0.4). Fig. [1] (left panel) 
shows the redshift distribution of 593 XMM-COSMOS 
AGN and the scaled random sample (~ 41000 random 
sources) which follows the red solid curve obtained by 
Gaussian smoothing. 

4. TWO-POINT STATISTICS 

A commonly used technique for measuring the spatial 
clustering of a class of objects is the two-point correlation 
function £(r), which measures the excess probability dP 
above a random distribution of finding an object in a vol- 
ume element d V at a distanc e r from another randomly 
chosen object (jPeebles Ill980h : 

dP = n[l + £(r)]dV (1) 

where n is the mean number density of objects. In par- 
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Fig. 3. — Projected AGN correlation function w p (r p ) computed 
at different r p scale (see label) as function of the integral radius 
iTmax- Horizontal lines show that the ACF saturates for irmax >40 
Mpc/h, which is also the minimum Umax at which w p (r p ) converges 
and returns the smaller error on the best-fit correlation parameter 
ro, with 7 fixed to 1.8. 

ticular, the auto-correlation function (ACF) measures 
the excess probability of finding two objects from the 
same sample in a given volume element. With a red- 
shift survey, we cannot directly measure £ (r) in physical 
space, because peculiar motions of galaxies distort the 
line-of-sight distances inferred from redshift. To separate 
the effects of redshift distortions, the spatial correlation 
function is measured in two dimensions r p and 7T, where 
r p and ir are the projected comoving separations be- 



tween the considered objects in the directions perpendic- 
ular and parallel, respectively, to the mean line-of-sight 
betwe en the two sources. Following iDavis fc Peeblesl 
(1983), r\ and r-i are the redshift positions of a pair of 
objects, s is the redshift-space separation (7*1 — ra), and 
/ = s(ri + r%) is the mean distance to the pair. The 
separations between the two considered objects across r p 



and 7r are defined as: 



si 



(2) 
(3) 



Redshift space distortions only affect the correla- 
tion function along the line of sight, so we esti- 
mate the so-called proj ected correlation function w p (r p ) 
(|Davis fc PeeblesllT983Tl : 



w p (r p ) 



£(r p ,n)diT 



(4) 



where £(r p ,7r) is the two-point correl ation function in 
term of r p and n, measured using the lLandv fc Szalavl 
(|1993l LS) estimator: 



f = \DD' - 2DR' + RR'} 

RR 1 



(5) 



DD', DR' and RR' are the normalized data-data, data- 
random and random-random number of pairs defined 
by: 



DD' = 



DR' 



RR' = 



DD{r p ,ir) 
nd{n d - 1) 

DR{t p ,-k) 
ridn r 

RR(r p , 7r) 
n r (n r — 1) 



(6) 
(7) 
(8) 



where DD, DR and RR are the number of data-data, 
data-random and random-random pairs at separation 
r p ± Ar p and it ± An and rid, n r are the total number 
of sources in the data and random sample, respectively. 
Fig. [1] (right panel) shows the number of pairs in redshift 
bins Az = 0.01 for the AGN sample. 

The LS estimator has been used to measure cor- 
relations in a number of surveys , for exam ple, SPSS 
(iZehavi et all 120051 : iLi et all 12006ft. DEEP2 (jCoil et al.1 
120071 120081) 7 AGESl Hickox et all l2009lh COSMOS 
( Gilli et al.ll2009D . If -K max = 00, then we average over 
all line-of-sight peculiar velocities, and w p (r p ) can be di- 
rectly related to £(r) for a power-law parameterization, 
by: 



w p (r p ) = r p ( — 
■ r p 



r yr(i/2)r[(7-i)/2] 



r(7/2) 



(9) 



In practice, we truncate the integral at a finite ir max 
value, to maximize the correlation signal. One should 
avoid values of ir max too large since they would add noise 
to the estimate of w p (r p ); if instead, n max is too small 
one would not recover all the signal. To determine the 
appropriate Tr max values for the XMM-COSMOS AGN 
correlation function, we estimated w p (r p ) for different 
values of ir max in the range 20-120 Mpc hr 1 . Besides, we 
determined the correlation length ro for this set of 7r ma2; 
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Fig. 4. — Left panel: Projected AGN ACF (black circles) compared to the auto-correlation of BL AGN (blue squares) and NL AGN 
(red triangles). The data points are fitted with a power-law model using the x 2 minimization technique; the errors are computed with a 
bootstrap resampling method. Right panel: The confidence contours of the power-law best-fit parameters ro and 7, for the whole AGN 
sample (black), for the BL AGN (blue) and NL AGN (red) sub-samples. The contours mark the 68.3% and 95.4% confidence levels 
(respectively corresponding to A\ 2 = 2.3 and 6.17) are plotted as continuous and dotted lines. 



values, by fitting w p (r p ) with a fixed 7=1.8 over r p in the 
range 0.5-40 Mpc hr 1 . In Fig. [3]we show the increase of 
the projected AGN auto-correlation w p (r p ) as a function 
of the integration radius ir max . The w p (r p ) values ap- 
pear to converge for ir max > 40 Mpc hr 1 . Therefore we 
adopt Tr max — 40 Mpc hr 1 in the following analysis, which 
is the minimum 7r mQX at which the correlation function 
converges. Such ir max selection returns the smallest error 
on the best-fit correlation parameter ro . 

5. PROJECTED AUTO-CORRELATION FUNCTION 

5.1. Standard Approach 

To estimate the AGN auto-correlation function £(r p , tt) 
using the LS formula (Eq. [5]), we created a grid with r p 
and 7r in the range 0.1-100 Mpc in logarithmic bins 
Alog(r p ,ir) — 0.2 and we projected £(r p ,7r) on r p using 
Eq. m 

In literature, several methods are adopted for error esti- 
mates in two-point statistics and no one has been proved 
to be the most precise. It is known that Poisson es- 
timators generally underestimate the variance because 
they do not account for the fact that the points are not 
statistically independent, i.e. the same objects appear 
in more than one pair. In this work we computed the 
errors on w„(r„) with a bootstrap resampling technique 
dCoil et alJl200aiHickox et alJl2009l iKrumpe etai]l2010t 
Capp elluti etaTlboiOf i. ~ 

The standard approach used to evaluate the power of 
the clu stering signal is t o fit w v (r v ) with a power-law 
model dCoil et all [2005 i Hickox et all [20091: iGilli et al.l 
l2009HKrumpe et alJl20ld : iCappelluti et al.ll2010h of the 
form given in Eq. [9[ using a x 2 minimization technique, 
with 7 and ro as free parameters. Fig. [4] (left panel, up- 
per quadrant) shows the projected AGN ACF, evaluated 
in the projected separation range r p = 0.5-40 Mpc hr 1 . 
The best-fit correlation length and slope and the corre- 



sponding la errors, are found to be r = 7.12±g;?| Mpc 
h- 1 and 7 = 1.81±g;g|. 

We estimated the projected correlation function of BL 
and NL AGN in the range r p = 0.5-40 Mpc hr 1 , 
as shown in FigU (left panel, lower quadrant). For BL 
AGN we found a correlation length r = 7.08t°;28 M P C 
h- 1 and 7 = 1.88±g;{$, while for NL AGN we measured 
r = 7-12±g;lg Mpc h' 1 and a flatter slope 7 = 1.69±g:g|. 
Fig. \Q(right panel) shows the power-law best-fit param- 
eters for the different AGN samples with the la and 2a 
confidence intervals for a two parameter fit, which corre- 
spond to X 2 = Xmin + 23 and X 2 = Xmin + 6 - 17 - 

We can estimate the AGN bias factor using the power- 
law best fit parameters: 



tpl 



os,AGn(z)/o-dm(z) 



(10) 



where as,AGN(z) is rms fluctuations of the density distri- 
bution over the sphere with a comoving radius of 8 Mpc 
hr 1 , auM^z) is the DM correlation function evaluated at 
8 Mpc h~ x , normalized to a value of aoM{z = 0) = 0.8. 
For a power-l aw correlation function this value can be 
calculated by ([Peebles Ill980h : 



(ff8,AGff) 2 = J r 2(7)( 



8Mpc/h 



(11) 



where J 2 (i) = 72/ [(3 - 7) (4 - 7) (6 - 7)2^]. As the lin- 
ear regime of the structure formation is verified only 
at large scales, the best-fit parameters ro and 7 are 
estimated fitting the projected correlation function on 
r p = 1 — 40 Mpc hr 1 . The la uncertainty of <T8,agn is 
computed from the ro vs. 7 confidence contour of the 
two-parameter fit corresponding to x 2 — Xmin + 2.3. 

5.2. Two-halo Term 

In the halo model approach, the two-point correlation 
function of AGN is the sum of two contributions: the 



6 



Allevato et al. 



TABLE 1 

Bias Factors and hosting DM halo masses 



(1) 


(2) 


(3) 


(4) 


(5) 


AGN 


< z > a 


bpL 




logM DM h 


Sample 




Eq. \10\ 


Eq. [TBI 




Total (593) 


1.22 


2-80t°;gg 


2.98 ±0.13 


13.23 ±0.06 


BL (354) 


1.55 


q 11+0.30 

— 1 22 


3.43 ±0.17 


13.14 ±0.07 


NL (239) 


0.74 


2.78l?;« 


2.70 ±0.22 


13.54 ±0.10 


X-unobs (184) 


1.12 


2.981°$ 


3.01 ±0.21 


13.33 ±0.08 


X-obs (218) 


1.30 


l-66±g-S 


1.80 ±0.15 


12.30 ±0.15 




Subsample at , 


5 < 1 




BL (70) 


0.57 


Z - L °-1.02 


2.32 ±0.26 


13.50 ±0.11 


NL (137) 


0.53 


1 fi«+ - 45 

l-o8_ 57 


1.40 ±0.15 


12.65 ±0.18 



a Median rcdshift of the sample. 
_ Typical DM halo masses based 
Ivan den Boschl 1 12003 1. 



ISheth et all ||2001D and 




Fig. 5. — Factor g as defined in Eq. 1211 estimated at the redshift 
of each AGN (black triangles). The data points are fitted by the 
function D\ (z) / Di (z = 0), where Z?i (z) is the growth function (see 
Eq. (10) in Eisenstein & Hu 1999 and references therein). The bias 
of each AGN is weighted by this factor according to the redshift z 
of the source. 

first term (1-halo term) is due to the correlation between 
objects in the same halo and the second term (2-halo 
term) arises because of the correlation between two dis- 
tinct halos: 



w A GN{r p ) = w\£ N {r p ) + w 2 AG h N (r p ) 



(12) 



As the 2-halo term dominates at large scales, we can con- 
sider this term to be in the regime of linear density fluc- 
tuations. In the linear regime, AGN are biased tracers 
of the DM distribution and the AGN bias factor defines 
the relation between the two-halo term of DM and AGN. 



W AGn( T p) ~ ^AGN w DM^ r p) 



2-ht 



(13) 



We first estimated the DM 2-halo term at the median 
redshift of the sample, using: 



dmM - 2 ^ 2 J 



P 2 -' l (k)k 2 



i(kr) 



kr 



,1k 



(14) 



i. io 



. ACF AGN 

b 4 wS i "(z=0) 



Fig. 



r p [Mpo h" 1 ] 

Projected AGN ACF (black circles) compared to 

~^(r v , z = 0) (dotted line), where the weighed bias b is de- 
fined in Eq. 1221 The shaded region shows the projected DM 2-halo 
term scaled by (b ± 8b) 2 . 

where P 2 ~ h (k) is the Fourier Transform of the linear 
power spectrum, assuming a power spect rum shape pa- 
rameter r = 0.2 and h = 0.7. Following lHamana et al.1 
(2002), we estimated Ci)M( r ) an d then the DM projected 
correlation (r p ) using: 



e D ^{r)rdr 



(15) 



Using this term, we can estimate the AGN bias simply 
dividing the projected AGN correlation function at large 
scale (r p > 1 Mpc h' 1 ) by the DM 2-halo term: 



h 2 

°AGN 



[waGnI - H r p)/wDM( r p)) 



a/2 



(16) 



and then averaging over the scales r p = 1 — 40 Mpc hr 1 . 
Table [TJ column 4 shows the AGN bias factors using this 
method, compared with the ones based on the power-law 
fits of the ACF (column 3) for the different AGN subsets. 
The two sets of bias values from the different approaches 
are consistent within la, but the errors on bpL are bigger 
consistently with the fact that the AGN ACF is not well 
described by a power-law. 

6. SOLVING FOR SAMPLE VARIANCE USING HOD 

The standard approac hes used in previous works on 
clust e ring of X-ray AGN dMullis et al Il2004t lYang et al ' 
2006 ] lG4Ui~eTaLl [20051: iCoil et all 120091: Iffickox et al 



2009; iKrumpe et alj|20l6t ICappelluti et al.ll2010l) to es 
timate the bias factors from the projected AGN ACF are 
based on the power-law fit parameters (method 1). This 
method assumes that the projected correlation function 
is well fitted by a power-law and the bias factors are 
derived from the best fit parameters rn and 7 of the clus- 
tering signal at large scale. 

Most of the authors dHickox et "all 120091 : IKrumpe et al.l 
[20101: ICappelluti et al.l l26TofT us ed an analytical expres- 
sion (as the one described in She th fc Tormenl 119991 : 
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Fig. 7. — Projected ACF of BL AGN (blue triangles, left panel) and NL AGN (red squares, right panel), compared to fe 2 w^Ji^p, z = 0) 
(dotted line), where the weighed bias b is denned in Eq. 1221 The shaded region shows the projected DM 2-halo term scaled by (b ± 5b) 2 . 
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Fig. 8. — Projected ACF of X-unobs AGN (darkblue open circles, left panel) and X-obs AGN (magenta diagonal crosses, right panel), 
0) 

2-halo term scaled by (b ± 5b)' 2 . 



compared to b 2 u>^, A J(rj,, z = 0) (dotted line), where the weighed bias b is defined in Eq. 1221 The shaded region shows the projected DM 



iSheth et al.l 120011 : iTinker et all 120051 ) to assign a char- 
acteristic DM halo mass to the hosting halos. The in- 
congruity of this approach is that the bias used is the 
average bias of a given sample at a given redshift. How- 
ever, the average bias is sensitive to the entirety of the 
mass distribution so that distributions with different av- 
erage masses, can give rise to the same average bias value. 
In the halo model approach the large scale amplitude sig- 
nal is due to the correlation between objects in distinct 
halos and the bias parameter defines the relation between 
the large scale clustering amplitude of the AGN ACF and 
the DM 2-halo term (method 2). 

In literature, the common model used for the AGN 
HOD is a three parameter model including a step func- 



tion for the HOD of central AGN and a truncated pow er- 
law satellite HOD (introduced bv lZehavi et al.ll2005l for 
galaxies). Here we assumed that all the AGN reside 
in central galaxies. This assumption is supported by 
iStarikova etaLl (|2010D . They found that X-ray AGN are 
predominantly located in the central galaxies of the host 
DM halos and tend to avoid satellite galaxies, fixing the 
limit to the fraction of AGN in non-central galaxies to 
be less than 10%. The same fraction of s at ellites galaxies 
hostin g AGN is suggested in lShenl ()2009f ). iShankar et"~aH 
(2010) modelled the measurements of quasar clustering 
derived in the SDSS (jShen et al.ll2009f ) and they verified 
that the predicted bias factors and the correlation func- 
tions are not altered including subhalos as quasar hosts. 
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TABLE 2 

Weighted Bias factors and hosting DM halo masses 



(1) 


(2) 


(3) 


(4) 


(5) 


AGN 


b 


z 


logMo 


&S0l a 


Sample 


Eq. [22] 


Eq. [231 







Total (593) 


1.91 


±0.13 


1.21 


13.10 


± 


0.06 


2.71 


± 


0.14 


BL (354) 


1.74 


±0.17 


1.53 


13.24 


± 


0.06 


3.68 


± 


0.27 


NL (239) 


1.80 


±0.22 


0.82 


13.01 


± 


0.08 


2.00 


± 


0.12 


X-unobs (184) 


1.95 


±0.21 


1.16 


13.30 


± 


0.10 


3.01 


± 


0.26 


X-obs (218) 


1.37 


±0.15 


1.02 


12.97 


± 


0.08 


2.23 


± 


0.13 




Subsample at 


z < 1 












BL (70) 


1.62 


±0.26 


0.63 


13.27 


± 


0.10 


1.95 


± 


0.17 


NL (137) 


1.56 


±0.15 


0.60 


12.97 


± 


0.07 


1.62 


± 


0.15 



a Bias estimated from M using lSheth et al.l J2001f) . 

A further consideration is that there is in practice no dis- 
tinction between central and satellite AGN in the 2-halo 
term that we used to estimate the AGN bias factor. 

We assumed a simple parametric form of the AGN halo 
occupation Na, described by a delta function: 



N A (M DM ) = f A 8(M_ 



DM 



M ) 



(17) 



where /a is the AGN duty cycle. It is clear that we are 
not considering the full HOD model, but we are assigning 
to all the AGN the same average mass of the hosting 
halos. The motivation is that X-ray AGN mainly reside 
in massive halos with a narrow distribution of the hosting 
halo masses. It's clear that this assumption is specific to 
AGN and e.g. is not applicable to galaxies. 

The AGN HOD descrived by ^- function is motivated by 
the results of lMivaji et al.1 (|2Qlll ) showing that the AGN 
HOD rapidly decrea ses at high halo masse s . In a ddition 
iMartini et all (|2009D and [Silverman et~aT1 (|2009ft found 
that AGN preferentially reside in galaxy groups rather 
than in clusters. 

The ^-function is the simplest possible assumption in the 
treatment of the sample variance, which is due to the 
variation in the amp litude of source counts dis tribution. 
It has been shown mlFaltcnbach er et al.l (|2010[ ). that the 
variation in the density field, which is responsible for the 
sample variance, can be replaced by the variation of the 
halo mass function. In terms of halo model, the bias 
factor as a function of the fluctuations A in the density 
field is expressed by: 



MA) 



J Mh N A (M h )b h (M h )n h (M h , A)dM h 
J M N A (M h )n h (M h ,A)dM h 



(18) 



where Na is the AGN HOD, b h (M DM ) is the halo bias 
and (.Mdm, A) is the halo mass function, which de- 
pends on the density field. On the other hand the sam- 
ple variance does not effect the AGN halo occupation. In 
Allevato et al. (in prep.) we confirm the assumption of 
constancy of the AGN HOD with the density field. 
When we assume that all AGN reside in DM halos with 
the same mass, Eq. [T5] becomes simpler: 



J Mh 5(M h ~ M a )b h (M h )n h (M h ,A)dM h 



f Mh S(M h - M )n h (M h , A)dM h 



b(M ) (19) 



The equation shows that when the AGN HOD is close 
to a (5-function, the variations in the density field only 



change the AGN number density and put more weight on 
AGN bias at the redshift of large scale structure (LSS), 
but do not change the bias of AGN inside th e structure. 
Our cl aim d iffers from the re sults presented in lGilli et al.1 
(|2005fl and (jGilli et al.ll2009h . They found that excluding 
sources located within a large-scale structures, the corre- 
lation length and then the bias factor strongly reduces. 
Such bias behaviour can be used to constrain more com- 
plicated shapes of the AGN HOD than a (5-function type 
distribution. 

However, even in the case of a (5-function HOD, we still 
need to consider the two effects which are often omitted 
in the clustering analysis: the LSS growth and the evolu- 
tion of the bias factor with z. Ignoring these effects can 
by itself lead to a difference in the results reported for 
the different AGN samples. 

The bias factor depends on the redshift as the struc- 
tures grow over time, associated with our use of a large 
redshift interval. For the i th source at redshift z%, we 
considered the bias factor corresponding to a halo mass 
M DM = M : 

h = b(M ,Zi) (20) 



whe re b(Mp,z) i s evalu ated using Ivan den Bosch! ([2002) 
and iSheth et al.l (|20Qll ). For each AGN at redshift z we 
estimated the factor g(z) defined as the square root of 
the projected DM 2-halo term at redshift z normalized 
to the projected DM 2-halo term evaluated at z = 0: 



9(*) = 



' WDM{z,r p ) 

w D m{z = 0,r p ) 



(21) 



averaged over the scales r p = 1 — 40 Mpc hr 1 . As the 
amplitude of the projected DM 2-halo term decreases 
with increasing redshift, g is a decreasing function of z 
(see fig. O, well described by the term D 1 (z)/ Di(z = 
0), where D\(z) i s the growth function (see eq. (10) in 
lEisenstein fc Hul (|2001[) and references therein). 
By accounting for the fact that the linear regime of the 
structure formation is verified only at large scales, we 
estimated the AGN bias considering only the pairs which 
contribute to the AGN clustering signal at r p = 1 — 40 
Mpc ft. -1 . We defined the weighted bias factor of the 
sample as: 



b{M ) 



(22) 



where bibj is the bias factor of the i and j source in 
the pair i—j, gigj is the g factor of the pair and N pa i r is 
the total number of pairs in the range r p = 1 — 40 Mpc 

h-K 

Similarly, we defined a weighted average redshift of the 
AGN sample, weighting the redshift of each pair for the 
g factor and the bias of the pair (bibj): 



_ Eij h ' b .,9'9.r 
z ■■ 



(23) 



where z v 



(zi + Zj)/2. Following this approach we 



can find the value of Mq that satisfies: 
h = b(M ) 
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where b\ is the square root of the projected AGN ACF 
normalized to the projected DM 2-halo term at z = 0: 



bi 



1 WAGN ( r p) 

wdm{z = 0,r p ) 



(24) 



averaged over the scale r p = 1 — 40 Mpc h -1 . 

By performing the test in narrow redshift intervals, we 
can study the dependency of the halo mass Mq on red- 
shift (see <JS]). Moreover with just a single measurement 
of the amplitude of the 2-halo term, one cannot constrain 
the AGN HOD. Already with several measurements sam- 
pling different density fields, the shape of the HOD can 
be linked to the LSS density-dependence of the bias. In 
addition, the 1-halo term of the AGN auto-correlation 
and AGN-groups cross-correlation can be used to dis- 
criminate between different HOD models, which will be 
argument of our following work. 

7. MEASUREMENTS 

The weighted bias factors b and redshifts z, and the 
corresponding DM halo masses Mq estimated for the dif- 
ferent AGN sub-samples using the method described in 
the previous section are shown in Table [2] 
Fig. HI [Hand E] show the ACF of the AGN, BL/NL AGN 
and X-unobs/obs AGN samples, compared to the term 



' J DM( r P' 



0) (dotted line), where the weighed bias 
b is defined in Eq. [22j The shaded region shows the 
projected DM 2-halo term scaled by (6± 5b) 2 . 

The AGN bias factor indicates that XMM-COSMOS 
AGN reside in halos with average mass logMo — 13.01 ± 
0.09[h~ 1 logMQ], characteristic of moderate-size poor 
groups, a result consistent with previous works on X- 
ray selected AGN that indicate that the typical DM halo 
mass hosting AGN is in the range 12.5 < XogMuM < - 
13.5[/i- x M ]. 

We found that BL and NL AGN which peak at z = 1.53 
and z = 0.82, present consistent bias factors which cor- 
respond to DM halo average masses logMg = 13.24 ± 



om[h- L M Q 

described in 



and 13.01 ± 1 08r fe~ 1 Af (3 ], respectively. As 
Brusa et al.l (|2010f ) . only a small fraction of 



the objects classified as NL AGN are located at z > 1, to 
be compared with 350 in the BL AGN sample. This is 
mostly due to the fact that high-redshift NL AGN are op- 
tically faint (typically / ~ 23 — 24) and have not been tar- 
geted yet with dedicated spectroscopic campaigns. Our 
results might be affected by the limitations in the ob- 
scured AGN classification, considering that some models 
on the evolution of the obscured AGN fra ction predict an 
incre ase of the fraction with the redshift (Hasingcr ct al. 
2008). In order to avoid the problem of different red- 
shift distribution in comparing BL/NL AGN cluster- 
ing amplitude, we selected for each sample a subset 
(BL AGN with 70 sources and NL AGN with 137) at 
Z ~ 0.6. At the same redshift we found that BL and NL 
AGN have a bias factor bsL = 1-62 ± 0.26 and b^L — 
1.56 ± 0.15, which correspond to average halo masses 
logM = 13.27 ± O.IO^M©] and 12.97 ± 0.07[/i" 1 Mo], 
respectively. 

Similar results have been obtained using X-unobs and X- 
obs AGN samples; unobscured AGN at z = 1.16 inhabit 
halos with average mass logMg = 13.30 ± O.lO[/i _1 M0] 
which is higher at 2.5 a level than the halo mass hosting 



obscured AGN (logM = 12.97 ± O.OS^Mq]), at sim- 
ilar redshift. 

In order to compare our results with previous works on 
the bias of X-ray selected AGN, we evaluated the bias 
factors corr e spond ing to the halo mass Mq at ~z using 
ISheth et al.l (|2001h as shown in Table col (5). 

Our results support the picture that at a given redshift, 
X-ray selected BL/X-unobs AGN reside in more massive 
halos compared to X-ray selected NL/X-obs AGN. This 
result would be expected if the two classes of AGN corre- 
spond to different phases of the AGN evolution sequence 



spor 
(Ho, 



Hopkins et"all[200l 120081 IHickox et al.ll2009f) . 

8. BIAS EVOLUTION AND CONSTANT MASS 
THRESHOLD 

In order to investigate the redshift evolution of the 
bias factor, we split the XMM-COSMOS AGN sample in 
three redshift bins. The sizes of the redshift bins have 
been determined such that there are more or less _the 
same number of objects in each bin. The values of 6, z 
and M for the total AGN sample are shown in Table 
03 The meaning of the table columns are: (1) sample; 
(2) number of sources; (3) bias parameter from the pro- 
jected DM 2-halo term, evaluated at the median < z > 
of the sam ple; (4) typical halo mass using Ivan den Bos"chl 
(p00l and ISheth etafl (pOOll ): (5) weighted bias of the 
sample; (6) weighted redshift of the sample; (7) Average 
DM halo mass; (8) Bias factor from Mo estimated using 
Sheth et al. (2001). 

We observed an increase of the AGN bias factor with 
redshift, from b(z = 0.92) = 1.80 ± 0.19 to b(z = 
1.94) = 2.63 ± 0.21 with a DM halo mass consistent 
with being constant at logM^n^M^] ~ 13.1 in each 
bin. These results support the picture that the bias of 
XMM-COSMOS AGN evolves with time according to a 
constant halo mass track at all redshifts z < 2. 

This conclusion, based on the analysis of the global 
XMM-COSMOS AGN sample, can however be affected 
by the fact that the relative proportions of BL and NL 
AGN are a strong function of redshift. In fact, since the 
XMM-COSMOS AGN sample is a flux limited sample, 
more luminous AGN are selected at high redshift and, 
also because of our magnitude limit, high- z sources in 
our sample are mainly BL AGN (see ^j. For this rea- 
son BL AGN sample could be analysed up to z ~ 2.25, 
while the maximum average redshift of the two redshift 
bins for NL AGN is z ~ 0.91. We found evidence of a 
strong increase of the BL AGN bias factor in four red- 
shift bins (see Table [3]), with a DM halo mass constant 
at logM Q [h~ l M Q ] ~ 13.28 at all redshifts z < 2.25. For 
NL AGN we estimated b(z = 0.62) = 1.59 ± 0.13 and 
b{z — 0.91) = 1.87±0.19, which correspond to a constant 
halo mass values logM^h -1 Mq] ~ 13.02. We split the 
X-unobs and X-obs AGN samples in two redshift bins 
up to z ~ 1.5 and we found that the bias of X-unobs 
AGN (X-obs AGN) evolves according to a constant halo 
mass consistent with the mass of BL AGN (NL AGN) 
hosting halos. Fig. \9\(left panel) shows the redshift evo- 
lution of the average DM halo mass Mg for all the AGN 
subsets. The horizontal lines represent the mean value 
of M for BL/X-unobs AGN (dashed-blue), NL/X-obs 
AGN (long dashed-red) and for the whole AGN sample 
(dotted-black). Fig. [9] (right panel) shows the redshift 
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Fig. 9. — Left Panel: DM halo mass M as a function of z for different AGN sub-samples (see legend). The horizontal lines show the 
mean value of Mo for BL/X-unobs AGN (dashed-blue), NL/X-obs AGN (long dashed-red) and for the whole AGN sample (dotted- black). 
Right Panel: Redshift evolution of the bias parameter fegoi °f different AGN sub-samples. The dashed lines show the expected b(z) of 
typical DM halo masses Mom based on[Shcth ct al. (2001]). The masses are given in logMjjM in units of Ii~ 1 Mq. BL/X-unobs AGN 
present a strong bias evolution with redshift with a constant DM halo mass logMo = 13.28 ± 0.07[/i — 1 Mq] up to z ~ 2.4. NL/X-obs AGN 
reside in less massive halos with logMo = 13.00 ± 0.06[/i — 1 Mq], constant at z < 1.5. 



TABLE 3 
Bias Evolution 



(1) 


(2) 


(3) 




(4) (5) 


(6) 


(7) 


(8) 


< z > a 


N 


b 2-h 




logM DM h b 


z 


logMo 


&soi c 






Eq. [16] 






Eq. US 


H^Mq 










All AGN 








0.80 


190 


2.70 ±0 


19 


13.48 ±0.10 1.80 ±0.19 


0.92 


13.12 ±0.06 


2.30 ±0.11 


1.30 


220 


3.10 ±0 


18 


13.21 ±0.10 2.14 ±0.18 


1.42 


13.07 ±0.08 


3.02 ±0.11 


2.07 


183 


5.18 ±0 


21 


13.30 ±0.11 2.63 ±0.21 


1.94 


13.18 ±0.08 


4.37 ±0.27 










BL AGN 








0.67 


70 


2.62 ±0 


20 


13.57 ±0.10 1.52 ±0.20 


0.70 


13.26 ±0.06 


2.16 ±0.25 


1.25 


108 


3.06 ± 


23 


13.24 ±0.08 2.02 ±0.23 


1.25 


13.21 ±0.08 


3.00 ±0.27 


1.71 


92 


5.37 ±0 


28 


13.60 ±0.08 3.57 ±0.28 


1.72 


13.32 ±0.08 


4.31 ±0.30 


2.46 


85 


6.82 ±0 


27 


13.41 ±0.10 4.02 ±0.27 


2.25 


13.28 ±0.10 


5.60 ± 0.42 










X-unobscured AGN 








0.65 


98 


2.46 ± 


17 


13.51 ±0.11 1.62 ±0.17 


0.80 


13.28 ±0.05 


2.34 ±0.18 


1.66 


86 


4.85 ±0 


18 


13.51 ±0.10 2.10 ±0.18 


1.54 


13.33 ±0.06 


3.90 ± 0.33 










NL AGN 








0.53 


137 


1.40 ± 


13 


12.65 ±0.12 1.59 ±0.13 


0.62 


13.01 ±0.05 


1.70 ±0.10 


1.02 


102 


2.11 ±0 


19 


12.88 ±0.15 1.87 ±0.19 


0.91 


13.04 ±0.07 


2.20 ±0.17 










X-obscured AGN 








0.73 


106 


1.80 ±0 


14 


13.01 ±0.11 1.51 ±0.14 


0.85 


13.03 ±0.06 


2.08 ±0.12 


1.84 


112 


3.51 ±0 


16 


12.94 ±0.13 1.96 ±0.16 


1.51 


12.95 ±0.06 


2.95 ±0.14 



a Median redshift of the sample. 

b Typical DM halo masses base d onfSheth et all (120011) and Ivan den Boscbl 12003) . 
c Bias estimated from A/q using Shcth ct al. (2001). 



evolution of the bias factors 6501 (Table [3J col (7)) for 
different AGN sub-samples. The dashed lines show the 
expected b(z) associated to the typical DM halo mass 
based on lSlieth et al.l (|2001l) . 

These results show that X-ray selected BL/X-unobs 
AGN reside in more massive DM halos compared to X- 
ray selected NL/X-obs AGN at all redshifts z at ~ 3cr 
level. This suggests that the AGN activity is a mass 
triggered phenomenon and that different AGN phases 
are associated with the DM halo mass, irrespective of 



redshift z. 

9. DISCUSSION 

9.1. Which DM halos host X-ray AGN? 

We have introduced a new method that uses the 2-halo 
term in estimating the AGN bias factor and that prop- 
erly accounts for the sample variance and the growth 
of the structures over time associated with our use of 
large redshift interval of the AGN sample. Using this ap- 
proach we have estimated an average mass of the XMM- 
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Fig. 10. — Bias parameter as a function of redshift for various X-ray selected AGN (black data points), X-ray selected BL/X-unobs AGN 
(blue data points) and X-ray selected NL/X-obs AGN (red data points) as estimated in previous studies and in this work according to the 
legend. Our results refer to the bias factor fosoi showed in Table[2]col (5). The dashed lines show the expected b(z) of typical DM halo 
masses Mom based on Sheth et al. (2001). The masses are given in logMnM m units of }i~ 1 Mq. 



COSMOS AGN hosting halos equal to logM [h,- 1 Mpc] = 
13.10 ± 0.06 which differs at ~ 1.6a level f r om th e typ- 
ical halo mass Mr, m base d on S heth et al.l (|2001l) using 
the methode 2 (see £|5.2[) . The difference between the 
standard method and our own method is also clear for 
the mass of BL and NL AGN hosting halos. We have 
found that BL AGN inhabit DM halos with average mass 
logM^h^Mpc] = 13.24 ± 0.06 at z = 1.53 while halos 
hosting NL AGN have average mass logM^h^ 1 Mpc] = 
13.01 ± 0.08. BL AGN reside in more massive halos 
than NL AGN also selecting two subsamples that peak 
at the same median redshift z ~ 0.6. We obtained simi- 
lar results using X-ray unobscured AGN at ~z~ — 1.16 and 
X-ray obscured AGN at z = 1.02 (logMolh' 1 Mpc] = 
13.30 ±0.10 and logMoih' 1 Mpc] = 12.97 ±0.08, respec- 
tively). 

Instead the typical halo mass based on iSheth et al.1 



(12001 using the AGN bias estimated with the method 2, 
strongly depends on the median redshift of the sample. 
According to the method 2, BL AGN at < z >= 1.55 
reside in less massive halos compared to NL AGN at 
< z >— 0.74, while the result is different selecting two 
samples of BL and NL AGN at the same < z >~ 0.5. 
Our results agrees with the majority of the recent studies 
of X-ray surveys which suggest a picture in which X-ray 
AGN are typically hosted in DM halos with mass in the 
range 12.5 < logM DM \h~ Y Mpc\ < 13.5, a t low (< 0.4) 
and high (~ 1) redshift. Star ikova et al.1 (|20100 found 
that Chandra /Bootes AGN are located at the center of 
DM halos with M > M min = 4 x 10 12 h" 1 M . This 
mass estimate represents a threshold value, since they 
are assuming a halo occupation described by a step func- 
tion (zero AGN per halo/subhalo below M m i n and one 
above it). Our approach, in terms of HOD, is completely 
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Fig. 11. — Bias pa r amete r as a functio n of redshift for op tically selected BL A GN from previous w orks fCroom ct al. 2005i), green-crosses; 
Porciani & Norbcrg (2006), green-stars; Shcn ct al. (2009), green-open squares; Ross ct al. (2009), gree-open triangle) and X-ray selected 
BL (blue triangles) and X-unobs (blue open-circles) AGN and NL (red squa res) and X-obs (re d crosses) AGN as estimated in this work. 
The dashed lines show the expected b(z) of typ i cal DM halo masses based on lSheth et al.| d2001|) and the dotted lines represent the passive 
evolution of the bias, as described in lFrv et al.l Ijl996l ). The bias of optically selected BL AGN evolves with redshift following an evolution 
at constant halo mass, with a typical mass which remains practically in the range logU DM ~ 12.5 - 13h- 1 M Q at all redshifts z < 2.25. 
X-ray selected BL/X-unobs AGN reside in more massive DM halos at all redshifts z < 2.25, according to a typical mass of the hosting 
halos constant over time in the range logMnM ~ 13 — 13.5/i — 1 Mq. The bias evolution of NL/X-obs AGN seems to indicate that they 
reside in DM halo mass logMjjM ~ 13/i _1 Mq constant at all z < 1.5. These results suggest the picture that X-ray selected BL AGN are 
triggered by secular processes as tidal disruption or disk instabilities instead of major mergers between gas-rich galaxies as confirmed by 
semi-analytic models and observations for optically selected quasars. 



different. We assume a halo occupation described by 5- 
function, supported by the fact that AGN only reside in 
massive halos (then the AGN HOD can be described by 
a narrow halo mass distribution at high mass values, but 
not by a step function). 

Fig. [10] shows the bias factors of X-ray selected AGN 
(black), BL/X-unobs AGN (blue) and NL/X-obs AGN 
(red) as estimated in different surveys (according to the 
legend). Our results refer to the bias factors 6501 showed 
in Tabled column (5). The dashed lines show the ex- 
pected b(z) ass uming a cons t ant ty pical DM halo mass 
M DM , based on lSheth etHI (|200ll i. 



The previ ous studie s of IGilli etHl ll2005h fo r 
the CDFN IGilli et al.l (|2009D . IMullis et all (|200l . 
lYang et al.l fl2006D for CLASXS AGN suggest the sce- 
nario in which the typical DM halo mass hosting X- 
ray selected AGN is lo gMnufo^Mm] ~ 13.5. The 
bias values measure d in |Gilli et all (12005!) o n CDFS, in 
iHickox et all (|2009f >. iCoil et al.l (|2009t l and lYang et al.l 
(2006) and in this work, correspond to a lower halo mass 
{logMoMlh -1 Mq] ~ 13). A possible explanation could 
be that at fixed redshift, the bias and then the mass of the 
hosting halo, depends on the luminosity of the sample. 
The same explanation might be applied to the results on 
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BL/X-unobs AGN. 

The bias estima t es at z < 1 for NL/X-obs AGN in 
ICappelluti et al.l ([2010) and in this work, seem to in- 
dicate that the mass of NL/X-obs AGN hosting halos is 
logM D M[h- x M Q ] ~ 13. 

9.2. Optically selected vs X-ray selected AGN 

We first found evidence of a redshift evolution of the 
bias factor of X-ray selected BL / X-unobs AGN (fig. [TT1 
blue data points) and NL/X-obs AGN (red data points). 
The bias evolves with redshift at constant average halo 
mass logMolh^Mol ~ 13.3 for BL/X-unobs AGN and 
logMolh^MQ] ~ 13 for NL/X-obs AGN at z < 2.25 
and z < 1.5, respectively. Fig. [TT] shows the expected 
b(z ) assuming a cons tant typical DM halo mass based 
on iSheth et al.1 (|200lD (dashed lin es) and the so called 
passive bias evolution (dotted lines iFrv et al1 [l996') . The 
observed bias evolution suggests an average halo mass 
of the hosting halos, constant over time in the range 
logM DAiih^ 1 Mq] = 13 — 13.5, instead of an evolution 
of the bias in a model in which objects are formed at a 
fixed time and their distribution evolves under the influ- 
ence of gravity. 

There have been several studies of the bias evolution of 
optical quasar with the redshift as shown in fig. [Til (green 
data point s), based on large survey samples such as 2QZ 
and SPSS (ICroom et al.ll2005l: IPorciani fc Norberdl2l)0a 
iShen et all 120091: iRoss et al.l I2009D . Since the quasar 
samples used in these clustering analysis are defined 
as spectroscopically identified quasars with at least one 
broad (FWHM>1000 km s _1 ) emission line, we refers to 
them as optically selected BL AGN. 

All the previous studies infer the picture that the 
quasar bias evolves with redshift following a constant 
mass evolution, with the average mass that can vary in 
the range logM r>M[h~ l Mq] ~ 12.5 — 13, may be depend- 
ing on the AGN sample luminosity as already suggested 
for X-ray selected AGN. The simplest interpretation ac- 
cording to the observed redshift evolution of the bias fac- 
tors is that 1) X-ray selected AGN whether BL/X-unobs 
or NL/X-obs AGN inhabit DM halos with mass higher 
than the mass of optically selected quasar hosting halos 
in the range z = 0.5—2.25; 2) X-ray selected BL/X-unobs 
AGN reside in more massive halos compared to NL/X- 
obs AGN for z — 0.6 — 1.6 and the discrepancy between 
the bias factors of the two samples increases with z; 3) 
the AGN activity is a mass triggered phenomena and the 
different AGN evolutionary phases are associated with 
just the DM halo mass, irrespective of the redshift z. 

9.3. External vs Internal Triggering 

The major merger of galaxies is one of the promis- 
ing mechanisms suggested to be responsible for fuelling 
quasars and in particular to be dominant for bright 
quasars at high redshift. Models of major mergers ap- 
pear to naturally produce many observed properties of 
quasars, as the quasar luminosity density, the shape 
and the evolution of the quasar luminosity function and 
th e large-scale quasar c lustering as a function of L and 
z (IHopkins et alJ [20081: IShenl [20091: IShankar et al.l I2009L 
2010; Shankar 20101: IBonoli et al. [ 12009ft . 
Clear evidence for higher incidence of mergers is seen 
among quasars ([Serber et alj|2006t IHopkins et £01 12006: 



— Shen (2009) at z=2 

Data Points: measured bias at 1.8 < z < 2.2 
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Fig. 12. — Predic t ed bia s as a function of luminosity, computed 
according to Shen (2009) fixing z = 2, compared to previous 
bias estimates at 1.8 < z < 2.2, for optically selected BL AGN 
and for XMM - COSM OS BL AGN. P oints are measuremen t s from 
IGroom et ail 120051 . green-crosses), Porciani & Norbcrg 120061 . 
green- star) , [Shen et al. (2009, gree n-ope n square). Ida Angela et af] 
120081 . green-circles ) , I Myers et al. (2007, green-squares) and our re- 
sult (blue triangle). For ease of comparison, all luminosities are 
converted to bolome tric luminosities using the corrections from 
Hopkins et al. (2007). The theoretical model which assumes a 
quasar phase triggered by major merger reproduces the results ob- 
tained for the bias of quasars, but can not reproduce the high bias 
factors found for X-ray selected BL AGN and then can not explain 
why optically selected quasars that have higher bolometric lumi- 
nosity compared to COSMOS X-ray selected BL/X-unobs AGN, 
reside in more massive halos. These differences suggest a switch 
to a different dominant mechanism for AGN triggering, from ma- 
jor mergers between gas-rich galaxies to secular processes as tidal 
disruptions or disk instabilities. 



iVeilleux et al .1120091) . Additionally a large fraction of lu- 
minous quasars at low redshift are associated with either 
morphologically disturbe d objects (jCanalizo fc Stockton! 
120011 : iGuvon et al.1 120061 ) . or early- type hosts with fine 
structure in their o ptical light distribution, indicative 
of pa st interactions ( Canalizo et al.| [2007: Bc nnert et al.1 
2008). In the local Universe, for instance, the study 
of the environm ent of Swift BAT Seyfert galaxies 
(|Koss et al.ll2010D appeared to show an apparent mergers 
~ 25% which suggests that AGN activity and merging 
are critically linked. Moreover it is believed that major 
mer ger dominates at high redshift and brigh t luminosi- 
ties ()Hasinger et ah 2008; Hopki n"s~et al.l[200 6). while mi- 
nor interaction or bar instabilities or minor tidal disrup- 
tions are important at low redshift (z < 1) and low lu- 
mino sities (Lbol ^5 10 44 erg s" 1 ) (Hop kins fc Henquisti 
[20091) . 

Our results on the bias evolution of X-ray selected 
BL/X-unobs AGN infer that these objects with Lbol ~ 
2 x I0 45 erg s^ 1 reside in massive DM halos Mum ~ 
2 x lO^Mm/i" 1 . Besides studies on BL AGN in the 
COSMOS field (jMerloni et al.l[20Tol : lTrTu"np et al.ll20H 
suggest that our sample is characterized by BH masses 
in the range Mbh = 10 7 — 10 9 M Q and Eddington ra- 
tio A > 0.01. Optically selected quasars from large sur- 
vey samples such as 2QZ and SDSS are high-luminosity 
quasars Lbol ^ 10 46 erg -1 with BH masses in the range 
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M B h = 10 8 - 10 10 M Q and A > 0.01. Clustering analysis 
of optical quasars have shown that they reside in DM 
halos with M DM ~ 10 12 M Q h~ x . 

Fig. [12] shows the predicted b ias as a func tion of lu- 
minosity computed according to IShenl ([2009) at z = 2. 
The theoretical model which assumes a quasar phase trig- 
gered by major mergers predicts an increasing bias with 
luminosity and reproduces the previ ous results obtained 
for optical quas ars at 1.8 < z < 2.2 (ICroom et~aLl (2005, 
green-crosses), iPorciani fc Norberel ([20061 green-star), 

IShen et all IpOOfll g reen-open square), Ida Angela et al.l 
(|2008l green-circles) JMvers et all ([2007L green-squares)). 
On the other hand the model can not reproduce the high 
bias factor found for X-ray selected COSMOS BL AGN 
(blue triangle) and then can not explain why optically 
selected quasars characterized by higher bolometric lu- 
minosity compared to X-ray selected COSMOS BL/X- 
unobs AGN, are found in less massive halos. These dif- 
ferences suggest a switch to a different dominant mech- 

a nism for AGN triggerin g. 

Hop kins fc Hernquist I ([2006D introduced a model for 
the fueling of low-luminosity AGN (Seyferts, with 
Lbol < 10 44 - 10 45 er 5 a" 1 and M BH < 1O 7 M ), 
which proposes AGN triggered by random accretion of 
gas via internal, secular processes. The stochastic accre- 
tion model and the merger-driven activity are fundamen- 
tally different, the former being determined by stochastic 
encounters with a cold gas supply in a quiescent system, 
the latter by the violent torquing of cold gas through- 
out entire galaxies into the galaxy center in major merg- 
ers. Accretion of cold gas in quiescent systems can ac- 
count for low luminosity Seyferts but can not explain the 
higher luminosities and the larger BH masses observed 
for XMM-COSMOS BL AGN. The high Eddington ra- 
tios at masses in the range M B h ~ 10 s — 10 9 M Q can not 
be maintained through this mode of accretion. 
Furthermore, this fueling mechanism predicts lower bias 
factors compared to the major merger picture for bright 
quasars, which is completely in disagreement with our 
results. 

Fueling by stellar winds or hot gas accretion may rep- 
resent yet a third qualitatively distinct mode of fueling. 
iCiotti fc Ostriker I (|1997L 120011 ) investigated the episodic 
AGN activity model in early-type galaxies, assuming 
at their center the presence of a massive BH growing 
with the accretion of matter and affecting the inflow 
through feedback. The duration of the single accretion 
event are extremely short but the maximum luminosi- 
ties reached during the accretion events can be of the 
order of Lbol ~ 10 46 — 10 47 erg a -1 , depending on the 
input parameters of the model. The central BH grows 
by episodic accretion up to a mass in the observed range 
(M ~ 10 8 - 5 - 1O 9 5 M ) in all giant ellipticals. 
On the other hand the observational consequence of this 
model is that the duty cycle is very low, typically of the 
order of 10~ 2 — 10 ~ 3 . This result implies a small fraction 
of giant ellipticals observed in an AGN phase, too low 
compared to the observed 10% of X-ray AGN residing in 
massive galaxies. 

In the AGN evolutionary model described in Hickox 
et al. (2009), optically bright quasars are hosted by on- 
going disk galaxy mergers and immediately precede an 
optically faint X-ray AGN phase, which evolves into an 



early-type galaxy. Following this evolutionary sequence, 
NL/X-obs AGN should be triggered in the first initial 
phase of vigorous star formation and obscured accretion 
which supports the scheme of NL AGN inhabiting halos 
with low typical masses logMr>M\h~ l M@\ ~ 12.5. An 
X-ray AGN phase immediately follows the quasar phase. 
Since DM halos grow and accumulate mass over time, X- 
ray AGN reside in more massive DM halos with typical 
mass logMpMih -1 Mq] ^13 — 13.5. This model predicts 
that X-ray AGN reside in more massive halos than QSO, 
but assumes a decline of the BH accretion rate from its 
peak in the quasar phase to M < 10 -2 M Edd or lower, 
which is in disagreement with the hi gh Eddington ratios 
found for XMM-CO SMOS BL AGN dMerloni et al.ll20ToT: 
iTrump et ai]|20ll . 

A plausible scenario requires that high-luminosity 
quasars (Lbol > 10 46 erg s _1 ) are triggered by external 
processes such as major mergers between gas-rich galax- 
ies with masses of the order of A/* ~ 10 10 A#q. Instead for 
BL AGN with Lbol ~ 2 x 10 i5 erg s" 1 , internal mech- 
anisms such as tidal disruptions or disk instabilities in 
massive galaxies (M* ~ lO n M0) might play a dominant 
role. 

The morphology of the AGN hosts galaxies provides 
an important clue into the mechanism that triggers their 
current AGN activity. It was observed that many AGN 
are not fueled by major mergers and only a small fraction 
of AGN are associated with m orphologically disturbed 
galaxies. Cistcrnas et al.l (|2010f ) analysed a sample of X- 
ray selected AGN host galaxies and a matched control 
sample of inactive galaxies in the COSMOS field. They 
found that mergers and interactions involving AGN hosts 
are not dominant and occur no more frequently than for 
inactive galaxies. Over 55% of the studied AGN sample 
which is characterized by Lbol ~ 10 45 er<? s~ 4 and by 
mass of the host galaxies M* > 1O 1O M0 are hosted by 
disk-dominated galaxies. This high disk fraction means 
that the lack of disturbed morphologies observed among 
the AGN hosts can not simply be due to a time lag be- 
tween merger activity and X-ray visibility and suggests 
that secular fueling mech anisms can be high efficie nt. 
It was also suggested bv iGeorgakaki s et al.l (|2009f ) that 
bar instabilities and minor interactions are more ef- 
ficient in producing luminous AGN at z < 1 and 
not only Seyfert galaxies and lo w -luminosity AGN as 
the [H opkins fc Hernquist I <|2006f >: [Hopkins & Hcnquist 
(2009) model p redicts. Besides several works on the AGN 
host galaxies ([Dunlop et al.l 120031: iGrogin et al.l 120051: 
Pierce et ail 120071 [Gabor et alj 120091 : iReichard et al.l 
2009t iTal et al.ll2009f) show that the morphologies of the 
AGN host galaxies do not present a preference for merg- 
ing systems. 

At the redshif t of o ur i nterest, recent findin gs of 
iSchawinski et al.l ([201 lh and iRosario et al.l (|2011| ). who 
examined a smaller sample of AGN at z ^ 2 in the ERS- 
II region of the GOODS-South field, inferred that late- 
type morphologies are prevalent among the AGN hosts. 
The role that major galaxy mergers play in triggering 
AGN activity at 1.5 < z < 2.5 was also studied in the 
CDF-S. Kocevski et al. (in prep.) found that X-ray se- 
lected AGN at z ~ 2 do not exhibit a significant excess 
of distorted morphologies while a large fraction reside in 
late-type galaxies. They also suggest that these late-type 
galaxies are fueled by the stochastic accretion of cold gas, 
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possibly triggered by a disk instability or minor interac- 
tion. 

We want to stress that our results by no means infer 
that mergers make no role in the AGN triggering. On the 
contrary, high luminosity AGN and probably a fraction 
of moderate luminosity AGN in our sample might be fu- 
elled by mergers. In fact, given the complexity of AGN 
triggering, a proper selection of an AGN sub-sample, us- 
ing for instance the luminosity, can help to test a partic- 
ular model boosting the fraction of AGN host galaxies 
associated with morphologically disturbed galaxies. 

Our work might extend the statement that for moder- 
ate luminosity X-ray selected BL AGN secular processes 
might play a much larger role than major mergers up to 
z ~ 2.2, compared to the previous z < 1, even during 
the epoch of peak merger-driven accretion. 

10. CONCLUSIONS 

We have studied the redshift evolution of the bias fac- 
tor of 593 XMM-COSMOS AGN with spectroscopic red- 
shifts z < 4, extracted from the 0.5-2 keV X-ray image 
of the 2deg 2 XMM-COSMOS field. We have described 
a new method to estimate the bias factor and the asso- 
ciated DM halo mass, which accounts for the growth of 
the structures over time and the sample variance. Key 
results can be summarized as follows: 

1. We estimated the AGN bias factor bsoi = 2.71 ± 
0.14 at z = 1.21 which corresponds to a mass of 
DM halos hosting AGN equal to logM [h~ 1 M Q ] = 
13.10±0.10. 

2. We split the AGN sample in broad optical emission 
lines AGN (BL) and AGN without optical broad 
emission lines (NL) and for each of them we con- 
sidered a subset with z = 0.6 and we found that BL 
and NL AGN present 6 SO i = 1.95±0.17 and & SO i = 
1.62 ± 0.15, which correspond to masses equal to 
logMoih^Mo] = 13.27 ± 0.10 and 12.97 ± 0.07, 
respectively. 

3. We selected in the hard band a sample of X- 
ray unobscured and X-ray obscured AGN accord- 
ing to the column density and we found that 
X-ray unobscured (X-ray obscured) AGN inhabit 
DM halos with the same mass compared to BL 
(NL) AGN with logMoih^Mol = 13.30 ± 0.10 
[logMolh^Mo] = 12.97 ±0.08). 



4. We found evidence of a redshift evolution of the 
bias factors for the different AGN subsets, corre- 
sponding to a constant DM halo mass threshold 
which differs for each sample. XMM-COSMOS 
AGN are hosted by DM halos with mass logM = 
13.12 ± O.O7[ft," 1 M ] constant at all z < 2, BL/X- 
ray unobscured AGN reside in halos with mass 
logMo = 13.28 ± O^f/i^M©] for z < 2.25 while 
XMM-COSMOS NL/X-ray obscured AGN inhabit 
less massive halos logM a = 13.00 ± 0.06[/i _1 M Q ], 
constant at all z < 1.5. 

5. The observed bias evolution for XMM-COSMOS 
BL and NL AGN at all z < 2.25, suggests that the 
AGN activity is a mass triggered phenomenon and 
that different AGN evolutionary phases are associ- 
ated with just the DM halo mass, irrespective of 
the redshift z. 

6. The bias evolution of X-ray selected BL/X-ray un- 
obscured AGN corresponds to halo masses in the 
range logMjjMih -1 M Q ] ~ 13 — 13.5 typical of 
poor galaxy groups at all redshifts. Optically se- 
lected BL AGN instead reside in lower density en- 
vironment with constant halo masses in the range 
logM DM [h- 1 M Q ] ~ 12.5-13 at all redshifts. This 
indicates that X-ray and optically selected AGN do 
not inhabit the same DM halos. 

7. The theoretical models which assume a quasar 
phase triggered by major mergers can not repro- 
duce the high bias factors and DM halo masses 
found for X-ray selected BL AGN up to z <~ 2.2. 
Our results might suggest the statement that for 
moderate luminosity X-ray selected BL AGN secu- 
lar processes such as tidal disruptions or disk insta- 
bilities play a much larger role than major mergers 
up to z ~ 2.2, compared to the previous z < 1. 
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